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ABSTRACT 

We present the Kepler light curve of KIC 4552982, the first ZZ Ceti (hydrogen-atmosphere pulsating 
white dwarf star) discovered in the Kepler field of view. Our data span more than 1.5 years with a 
86% duty cycle, making it the longest pseudo-continuous light curve ever recorded for a ZZ Ceti. This 
extensive data set provides the most complete coverage to-date of amplitude and frequency variations 
in a cool ZZ Ceti. We detect 20 independent frequencies of variability in the data that we compare 
with asteroseismic models to demonstrate that this star has a mass > 0.6 Mq. We identify a 
rotationally split pulsation mode and derive a probable rotation period for this star of 17.47 ± 0.04 
hr. In addition to pulsation signatures, the Kepler light curve exhibits sporadic, energetic outbursts 
that increase the star’s relative flux by 2-17%, last 4-25 hours, and recur on an average timescale of 
2.7 days. These are the first detections of a new dynamic white dwarf phenomenon that we believe 
may be related to the pulsations of this relatively cool (Teff = 10, 860 ±120 K) ZZ Ceti star near the 
red edge of the instability strip. 

Subjeet headings: white dwarfs, stars: oscillations, stars: activity, stars: individual (KIC 4552982, 
WD J191643.83±393849.7) 


1. INTRODUCTION 

As they cool, white dwarfs — the endpoints of more 
than 97% of stars in our Galaxy — evolve through in¬ 
stability strips on the H-R Diagram where they pulsate 
due to co nvective driving (|Brickhill||1991| |Goldreich & 
Wu 1999). ZZ Ceti variables have partially ionized hy¬ 
drogen atmospheres and pulsate in the surface temper¬ 
ature range 12,600 K > Teff > 10,800 K near the mean 
hydrogen-atmosphere (DA) white dwarf mass of 0.6 Mq 
( Tremblay et al.|2013[ ). The characteristics of the pulsa- 
tion modes that are excited in stars are determined by 
the specifics of their internal structures. The tools of 
asteroseismology enable us to interpret measurements of 
white dwarf brightness variations to potentially constrain 
their masses, radii, compositions, chemical stratification, 
equations of st ate, rotation, crystallized fractions, etc, 
(see reviews bylWinget fc Keple r|20'Q8l lEontaine fc Bras- 
sard||2008t |Althaus et al.||2010|). 'I'his method allows us 
to conduct important investigations into the behavior of 
matter under the extreme physical conditions of white 
dwarf interiors that are beyond what is accessible for 
study in terrestrial laboratories. 

Asteroseismology of white dwarf stars is conducted 
primarily through Eourier analysis of photometric light 
curves. While ground-based observations of pulsating 
white dwarfs have been collected since 1964 (jLandolt 
1968), the terrestrial vantage point comes with its disad- 


keatonb@astro.as.utexas.edu 

^ Department of Astronomy, University of Texas at Austin, 
Austin, TX-78712, USA 

^ McDonald Observatory, Fort Davis, TX-79734, USA 
^Department of Physics, University of Warwick, Coventry- 
CV4 7AL, United Kingdom 

Penn State Worthington Scranton, Dunmore, PA-18512, 
USA 

^ Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 
B-3001 Heverlee, Belgium 


vantages. Aliasing caused by daily gaps in the data and 
breaks due to inclement weather can introduce artifacts 
to Eourier transforms (ETs). Since spectral resolution 
and signal-to-noise improve with a longer observational 
baseline and more complete coverage, observations span¬ 
ning multiple nights are often required to resolve the in¬ 
dividual frequencies of a rich pulsation spectrum. Efforts 
to surpass these limitations include extended, global ob¬ 
servations with networks of telescopes distributed in lon¬ 
gitude, with the greatest contributions to white dwarf 
science coming from the W hole Earth Telescope collab¬ 
oration ( Nather et al.||199Q ). 

More recently, the Kepler space mission has enabled 
asteroseismic research of unprecedented quality by ob¬ 
taining extended time series photometr y of a consistent 
field with a high duty cycle (see, e.g., Gilliland et al. 


2010a Christensen-Dalsgaard & Thompson' 201 ip . Be¬ 


sides the presently described work, detailed asteroseismic 
studies of pulsating white dwarfs enabled by the origi¬ 
nal Kepler mission include analysis of a V777 Her star 


Ostensen et al.| 20111 

i Bischoff-Kim & Ostensen 

2011 

Corsico et al. 2012 I 

lischoft-Kim et al. 2014 

P and an- 

other ZZ (Jeti (KiC 11911480; Greiss et al. ! 

^014). In 


its present two-wheel configuration, Kepler has also pro- 
vided asteroseismic d ata on the ZZ Ceti variable CD 1212 
(Hermes et al. 201 4|) and a ZZ Ceti w ith an M dwarf 
binary companion (jHermes et al.||2015). Kepler will ob¬ 
serve addition al white dwarf pnls ators in upcoming K2 
m ission fields dHo w ell et al.||2Q14P. 

Hermes et ah' j2Ull| ught and discovered the 


first ZZ Ceti in the original Kepler field of view: 
WD J191643.83±393849.7. They gathered ^21 hr of 
time-series photometry on the 2.1-meter Otto Struve 
telescope at McDonald Observatory and identified seven 
frequencies of brightness variability, though with admit- 
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Fig. 1.— Representative segments of the Kepler light curve of KIC 4552982 as a function of days since the start of Q14 observations. 
The top panel shows the full Q14 light curve. The one-month shaded region in the top panel is expanded in the middle panel. The 
one-week shaded region in the middle panel is expanded in the bottom panel. The solid line is the light curve smoothed with a 3-hour-wide 
Epanechnikov (inverted parabola) kernel. The point-to-point scatter dominates the pulsation amplitudes in the light curve, so pulsations 
are not apparent to the eye. The white dotted line marks the significance criterion for our outburst detection algorithm, and outbursts 
determined to be significant are highlighted. We discuss this algorithm and these outburst events in Section 


tedly large uncertainties. They also obtained four low- 
to medium-resolution spectra for the white dwarf and fit 
the Balmer line profiles to models to determine its val¬ 
ues of Teff = 11,129 ± 115 K, \ogg = 8.34 ± 0.06, and 
= 0.82 ± 0.04 Mq. Now equipped with corrective 
terms that take into account the effects of 3-dimensional 


convection (Tremblay et al. 2013), we revise the spec- 
troscopically derived values to 4eff = 10,860 ± 120 K, 
log^ = 8.16 ± 0.06 in thi s work, and interpol ate the 
model cooling sequences of Renedo et al. (2010) to get 
= 0.69 ± 0.04 Mq. These properties place the white 
dwarf at the e mpirical red (cool) edg e of the ZZ Ceti 
instability strip ( Tremblay et al.|[2013 ). 

This target at 19^ 16^43®.83, 39°38'49.7" was assigned 
Kepler ID KIC 4552982 and was observed by the Kepler 
spacecraft at short cadence from Qll until the second 
reaction wheel failure during Q17. The resulting data 
provide the longest (^20-month) pseudo-continuous light 
curve of a ZZ Ceti ever obtained. Besides resolving a 
rich pulsation spectrum, the Kepler light curve revealed 
a surprising outburst phenomenon unlike any previously 
studied white dwarf behavior. 

In Sectionwe describe the methods used to optimally 
reduce the Kepler light curves. In Section we charac¬ 
terize and discuss the nature of the energetic outbursts 
recorded for the first time in the Kepler photometry. Our 


asteroseismic analysis makes up Sectionj^ where we mea¬ 
sure the pulsational properties of the star to constrain 
the stellar mass, hydrogen layer mass and stellar rota¬ 
tion rate. We summarize our findings and conclude in 
Section [H 

2. Kepler PHOTOMETRY AND DATA REDUCTION 
Eollowing the discovery by [Hermes et al.' (2011) that 


WD J191643.83+393849.7 (KlC 4552982) pulsaTes, this 
target was prioritized by the Kepler Asteroseismic Sci¬ 
ence Consortium (KASC) for short-cadence m onitoring 


by the Kepl er spacecraft (58.85 s image“^; Gilliland 
et al.||2010b|). It has a recorded Kepler magnitude of 
Kp = 17.9. The object was observed from Quarter 11 
(Qll) through Q17 (2011 September 29 to 2013 May 11). 

The Kepler data were reduced with the PyKE soft¬ 
ware package ( Still fc Barclay|20 12) following the method 
described by Kinemuchi et al.| (|2012 ). We defined cus¬ 
tom apertures for light-curve extraction from the Target 
Pixel Piles to maximize target signal-to-noise. We then 
masked out the outburst events (described in Section 3) 
from the raw light curves and found that linear least- 
squares fits of the top six Cotrending Basis Vectors from 
each quarter characterized and allowed us to correct for 
the systematic instrumental trends in the light curve. Ei- 
nally, we excluded points that were flagged for question- 
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able data quality by the Kepler pipeline, and we man¬ 
ually removed statistical outliers falling > 5cr from the 
local, 10-d median, ensuring that no clipped points were 
related to astrophysical photometric variations. Our fi¬ 
nal light curve contains 746,916 data points and has a 
86% duty cycle. The representative Q14 light curve that 
we acquired after normalizing and combining each month 
of data is provided in Figure 

3. OUTBURSTS IN THE Kepler LIGHT CURVE 

With a mean standard deviation from measurement 
noise of 1.8% dominating the pulsation amplitudes in 
the Kepler light curve, its visually striking features are 
the occasional large brightness enhancements that occur 
throughout the 20 months of observations (visible in Eig- 
ure[^. We see no evidence of a faint, red companion, and 
argue these observations likely mark the first detection 
of a new white dwarf outburst phenomenon. 


3.1. Outburst Charaeteristies 

We define an automatic detection algorithm to locate 
significant outbursts in the light curve. Each month of 
data is smoothed with a 3-hour-wide Epanechnikov (in¬ 
verted parabola; Epanechnikov 1969) kernel to reduce 
scatter while retaining the outburst signatures. Out¬ 
bursts are identified where the light curve exceeds a sig¬ 
nificance threshold of 4x the standard deviation of the 
smoothed flux for 30 consecutive minutes. We interpret 
the target to be in outburst until the smoothed flux re¬ 
turns to the median value of 1.0. We iteratively search 
for outbursts by recalculating the significance threshold 
from “quiescent” (outside known outbursts) portions of 
the light curve until no additional outbursts are identi¬ 
fied. The values that we adopt for this detection scheme 
are tuned such that the algorithm identifies all outbursts 
that are obvious to the eye without selecting regions that 
do not stand up to human scrutiny. Our search yields 
178 events of significant (peak 2 — 17%) brightness in¬ 
creases that typically last 4 — 25 hours. We display 
the histogram of delay times between detected outbursts 
in Eigure These outbursts appear to occur stochas¬ 
tically in time with an average delay of 2.7 days and a 
standard deviation of 1.3 days. 



Delay Between Outbursts (days) 


Fig. 2.— Histogram of delay times between successive outbursts 
that meet our detection criteria. The mean recurrence timescale is 
2.7 days and the standard deviation is 1.3 days. 
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Fig. 3.— Autocorrelation of the Q11-Q17 Kepler light curve of 
KIC 4552982. The positive autocorrelation coefficients at short 
(< 8 hour) time lag simply demonstrate that outbursts often per¬ 
sist to this duration. The range of time lags that produce negative 
autocorrelation coefficients is interesting, as it shows that an out¬ 
burst is less likely to happen within 2 days of another. Beyond 
this recharge time, outbursts are consistent with being random 
(the solid lines are the 95% confidence intervals for random data). 


These timescale properties are reflected in the auto¬ 
correlation of the light curve, shown in Eigure The 
autocorrelation function depicts how likely it is to mea¬ 
sure excesses in flux separated by different lags of time. 
The light curve was smoothed in 30-minute bins to av¬ 
erage over any correlation from pulsations. The solid 
horizontal lines mark the 95% confidence thresholds of 


±2/yWe , where is the effective sample size (Chatfield 
2004). We approximate the effective sample size from 


the number of light curve bins compared at each time 
lag, multiplied by (1—ACEi)/(!-!-ACEi), where ACEi is 
the autocorrelation coefficient at the smallest time lag. 
Eor a random time series, 95% of autocorrelation coef¬ 
ficients would be expected to fall within these bounds. 
We see positive autocorrelation at < 8-hour time lag as 
this is the characteristic duration of an outburst. Beyond 
an apparent recharge time of roughly 2 days, outbursts 
are consistent with occurring at random. The outbursts 
appear to be aperiodic, so we cannot formulate an event 
ephemeris to predict the timing of future outburst events. 

Eigure [^displays regions of the light curve surrounding 
a few representative events in more detail. The top panel 
shows the most energetic event found in the light curve, 
the second panel shows a median-energy event, the third 
panel shows a multi-peaked event, and the bottom panel 
displays the mean event profile that is representative of 
these events as a whole. We calculated this mean profile 
by aligning the moments of peak brightness of all the 
detected events and averaging points into 58.85-s time 
bins. 

To quantify the energies of these outbursts in the Ae- 
pler bandpass, w e measure their equivalent durations 
(Gershberg 1972). These are measured in the same 
way as spectral line equivalent widths: we integrate the 
flux above the local mean for the duration of the event, 
normalized to the mean flux level. Assuming that the 
outbursts originate from the target star, this results in 
a value with units of time equaling the duration that 
the star would shine in quiescence to output the same 
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Fig. 4. — Examples of outburst events. Top: most energetic. Sec¬ 
ond: median energy. Third: multi-peaked. Bottom: mean profile. 
The black curves are smoothed with a 3-hour-wide Epanechnikov 
kernel. The regions determined as belonging to the outburst by 
our algorithm are indicated with darker gray points. 



Equivalent Duration (minutes) 


Fig. 5.— Histogram of equivalent durations for the 167 outbursts 
observed in their entirety. 


amount of energy in the observed bandpass as the mea¬ 
sured excess from the brightening event. A histogram of 
the measured equivalent durations of the 167 outbursts 
that were recorded without interruption from gaps in the 
data is displayed in Figure and the continua used for 


the example outbursts are included in Figure The 
median equivalent duration of detected outbursts is 5.5 
minutes (the corresponding outburst is displayed in the 
second panel of Figure [4). Our detection algorithm is 
tuned to detect large outbursts and is likely incomplete 
in identifying the lowest-energy events. 


3.2. Discussion of Outbursts 

The outburst characteristics that we measure are un¬ 
like any previously studied white dwarf behavior. We 
argue here that these events originate from the white 
dwarf target KIC 4552982 and mark the first detection 
of a new astrophysical phenomenon. The data in hand 
are insufficient for us to identify the physical mechanism 
that drives these events. While we suggest possible con¬ 
nections to the evolutionary state of KIC 4552982, more 
intensive observational and theoretical explorations into 
their nature are left for future work. 

Each quarter of Kepler data was separated by rolls of 
the spacecraft that positioned the KIC 4552982 stellar 
image on a repeating sequence of four different CCD de¬ 
tectors. Since the outbursts are prevalent in all quarters 
of Kepler data, they cannot be instrumental artifacts of 
any one of the detectors. Furthermore, these are not 
widespread systematic artifacts since they do not cor¬ 
relate with any of the Kep/er-provided Cotrending Basis 
Vectors that catalog trends common amongst light curves 
for many objects. No other short-cadence targets show 
this phenomenon. We also find no correlation between 
the outbursts and the short-cadence light curve that we 
extracted for KIC 4552992 - a star with Kp = 18.934 
mag that is separated from KIC 4552982 by 12.5” and 
fell mostly within the same CCD subregions read out for 
the target. We note that this other source is well sepa¬ 
rated from our target and did not contribute flux to our 
custom photometric apertures. _ 

From the discovery paper (Hermes et al.|2011 ), we have 
five low- to medium-resolution spectra of KiC 4552982 
that cover a combined wavelength range from 4500 to 
7200 A. None show evidence of a companion from ei¬ 
ther radial velocity variations or spectral features that 
are not common of a DA (hydrogen atmosphere) white 
dwarf in this wavelength regime. Since the discovery pa¬ 
per on KIC 4552982, new 3-dimensional convection sim¬ 
ulations have mapped the parameters determined from 
1-dimensional spectroscopic model fits to values that 
largely res olve a mass discrepan cy in the 1-dimensional 
approach (Tremblay et al. We apply these correc¬ 

tive ter ms and revise the spe ctroscopically derived val¬ 


ues from Hermes et al. (2011) to Teff = 10,860 ± 120 K, 
log^ = 8.16 ± 0.06. The white dwarf radiu s, as inter¬ 
polated fr om the model cooling sequences of |Renedo et| 
al. (20101), enables the conversion to = 0.69 ± 0.04 
Mq. These parameters maintain that KIC 4552982 
is one of the coolest known ZZ Ceti pulsators with a 
mass near or slightly above th e peak of the DA w hite 
dwarf mass distribution (e.g.. Falcon et al. 2010 find 
(Mda) = 0-647lo.oi4 gravitational redshift 

measurements). 

Apparent photometric magnitudes for KIC 4 552928 
are availabl e from the Kepler-INT Survey (KIS; [Greiss 


et al. 2014) in U^g^r, and i filters. We include these 


magnitudes in Table along with the J-band magni- 
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TABLE 1 

Apparent photometric magnitudes 
OF KIC 4552982 


Filter 

Vega magnitude 

AB magnitude 

U 

17.362 


0.007 

18.15 


0.01 

9 

17.755 


0.005 

17.68 


0.01 

r 

17.677 


0.007 

17.84 


0.01 

i 

17.565 


0.009 

17.94 


0.01 

J 

17.76 


0.03 

18.70 


0.03 


TT-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-[“ 


I 

(D 


Koester White Dwarf Model 
lllOOK, log^ = 8.25 


iS 


I li t f\ 

'! r\ 

Ij !■ \ 


\ 

A 


_L 


_L 


_L 


_i_I_I_i_ 


4000 6000 8000 10000 

Wavelength (A) 


12000 


Fig. 6. — The spectral energy distribution of KIC 4552982 in 
t/, 5 f, r, z from the KIS survey and a J-band magnitude from the 
UKIRT J-band Public Archive. The synthetic spectrum IKoester] 
(|2010|) near the spectroscopically determined atmospheric param- 
eters is plotted for reference. We see no infrared excess, allowing 
us to rule out a cool dwarf star companion to the white dwarf as 
the source of the observed outbursts. 

tude acquired from the UKIRT public archiv^ All of 
these are reported in the Vega photometric system. To 
compare these magnitudes in a spectral energy distribu¬ 
tion (SED), we convert from the Vega system to physical 
A B magnitudes using the c onversion factors calculated 
by [Blanton fc Ro weis (|2QQ7|) for the U, g, r, i filters, and 
the conversion of Hewett et al. (2006) for the J filter. 
We plot the SED of KiC 4552982 against a synthetic 
spectrum of a Teff = 11,000 K, log^ = 8.25 white dwarf 
(Koester 2010) for reference in Eigure[^ We note that 
the photometry appears consistent with me model out to 
the infrared, demonstrating the absence of any cool main 
sequence companion that we might entertain as the pos¬ 
sible source of the brightening events in the Kepler data. 
We did not deredden the magnitudes, which might ac¬ 
count for the slight discrepancy in the U-band. 

It is not likely that the outbursts are transient events 
from some line-of-sight flare star, rather than intrinsic 
to the pulsating white dwarf. While the slightly shorter 
rise time than decay time of the mean outburst profile 
in the bottom panel of Eigure is qualitatively remi¬ 
niscent of flares observed from low-mass main sequence 
stars, the observed flux increases are not nearly as rapid 


as th e impulsive phases of c lassical stellar flares (Moffett 
19741 iBenz & Gfidel||2010|). Extensive studies of hares 


from early-type M dwarfs (M 4 and earlier; Hawley et al. 


2014 Davenport et al.|2Q14 ) in the short cadence Aep/er 


data and fro m later-type M dwarfs from ground-based 
photometry (|Hilton||2011|) do not find such frequent, en¬ 
ergetic flares as we detect in our data. The durations 
and pulse shapes of flares from K dwarfs in the Kepler 
long-cadence data also do not match our observations 


(Walkowicz et al. 2011). That the quiescent flux of a 


line-ol-sight hare star m the Kepler bandpass must be 
orders of magnitude less than the white dwarf target to 
not be detected in the observed spectra or SED means 
that the already improbable energies/peaks/durations of 
the outbursts are only extreme lower-limits if they do not 
originate from the white dwarf. 

To constrain whether the outbursts could be coming 
from any spatially offset faint transient source — flare 
star or otherwise — we searched for correlations be¬ 
tween the photometric centroids within the Kepler im¬ 
ages (plate scale ~ 4" pixel” a nd the occurrences o f 


(2013) 


lotome- 


outbursts using the approach of [Bryson et al. 

Our analysis was inconclusive. Although the p 
try shows significant centroid offsets when comparing the 
light curve during outbursts to quiescent moments, the 
shifts are in inconsistent directions on the sky in different 
quarters. When comparing just the higher-flux to lower- 
flux frames from quiescent segments of the light curve, an 
identical centroid analysis again shows significant shifts 
of comparable magnitude in identical directions. The 
measured centroid shifts of as much as 0.04" are likely 
dominated by variations in pixel sensitivity rather than 
the presence of any transient offset by more than of order 
0 . 1 ". 

Some white dwarfs host surface magnetic fiel ds locally 
as strong as ^ 1000 MG (e.g., PG 1031-1-234; Latter et 


al. 


ca. 


1987). We do not undertake any attempt to physi- 


y model the output of magnetic reconnection flares 
in white dwarfs for this analysis; however, the available 
spectra of KIG 4552982 do not show line splitting that 
would evidence that it is highly magnetic ( resolved Zee- 
man splitting is expect ed for H > 10^ Gauss; [Wickrama^ 
inghe fc Eerrario||2000 ). Since this first object to exhibit 


this behavior is not nearly one of the more strongly mag¬ 
netic white dwarfs studied, it is difficult to ascribe the 
outbursts to magnetic processes. Since the timescales of 
dynamic events like flares should scale roughly with the 
dynamical timescale (free-fall time), we do not expect to 
see flares of multi-hour duration on a white dwarf with 
a ^ 1 s dynamical timescale. 

We collected 42 hours of follow-up, higher signal-to 
noise photomet ry with the Argos photometer (Mukadam 
& Nather|2005|) on the McDonald Observatory 2.1-meter 
Otto Struve telescope over eight nights through a HG40 
filter in 2013 August with hopes of catching the system 
in outburst. Unfortunately we made no such detection. 
Given the distribution of observing windows when we 
were able to get useful data, there was a 47% chance 
that we would have caught the peak of an outburst. We 
also do not detect this phe nomenon when reana lyzing the 
o riginal discovery data o f 


De Marco et al. 


Hermes et al. 


( 2011 ) 


(2015) analyze the Kepler 


http://keplerscience.arc.nasa.gov/ToolsUKIRT.shtml 


ight curve 

of the central star of planetary nebula Kn 61, which ex¬ 
hibits brightening events every 2-12 days, each lasting 
1-2 days with amplitudes of 8-14%. These events are 
similar to the outbursts on KIG 4552982, at least in that 
they are also unlike previously studied behavior. While 
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Fig. 7.— Top: the thermal timescale as a function of depth in 
our model. The 10-hr average outburst duration timescale and the 
2.7-day recurrence timescale are marked with vertical dotted lines. 
Bottom: the fractional flux carried by convection through the static 
model atmosphere. The relevant outburst timescales are of order 
the thermal timescale at the base of the convection zone, which 
is the timescale suspected to be most relevant to the cessation of 
pulsations at the cool edge of the instability strip. 


the authors are unable to conclusively demonstrate a 
mechanism for observed outbursts in Kn 61, they ex¬ 
plore accretion as a possible energy source. Extending 
this analysis, they equate the energy of a median out¬ 
burst on KIC 4552982 to 5 x 10“^^ Mq of accreted ma¬ 
terial. We note with some interest that this is within 
the range of measured asteroid masses (Hilton 2002), 
but there is currently no known dynamical mechanism 
to drop circumstellar debris onto a white dwarf with this 
fairly regul ar and relatively rapid recurrence timescale 
( Jura|[MQ8 ). Assuming a total asteroid population mass 
of order the Solar Sy stem asteroid belt mass (1.8 x 10^^ 
g;|Binzel et al.||2000 ), this rapid accretion would deplete 
the mass reservoir m a mere ^1000 yr. With an out¬ 
burst frequency far shorter than the calcium diffusion 
timescale of ^40 yr for a Teff = 10,8 60 K, log^ = 8.16 


white dwarf (Koester & Wilken 


2006), we would expect 


such rapid accretion to cause strong absorption at the 
calcium H & K doublet. We do not detect th is signature 
in the available spectra (jHermes et al.||2011|). 

With Teff = 10,860 A 120 K, Kllrte^82 is o ne of 
the coolest ZZ Get is known ([Tremblay et al.||2013|). The 
mechanism that shuts down pulsations at the empirical 
cool edge of the ZZ Ceti instability strip is not fully un¬ 
derstood, but may be related to the thermal timescale 
at the base of t he convection zone exceed ing some crit¬ 
ical value (e.g.. Van Grootel et al. 2013). To explore 
this timescale in KIC 455298 2, we use the Wars aw-New 
Jersey stellar envelope code (Pamyatnykh 1999) to cal¬ 
culate a static pure-hydrogen atmosphere model that 
matches the measured spectroscopic parameters with the 
ML2/a = 1.0 mixing-length prescription for convection. 


We display the run of the thermal timescale and the con¬ 
vection profile in the outer 10“^ of the star by mass in 
Figure We note that the thermal timescale at the 
base of the convection zone is of the same order as the 
outburst duration and recurrence timescales (and much 
longer than the ^1 second dynamical timescale). If this 
is more than just a numerical coincidence, then it is possi¬ 
ble that this outburst behavior is a more general property 
of ZZ Ceti variables as they evolve out of the instability 
strip. 

While more than 180 ZZ Cetis have been previously 
discovered and studied, there has been a slight bias to¬ 
wards the detection of pulsators near the hot edge of the 
instability strip. This is because their pulsations are ex- 


lutionary cooling (e.g., Kepler et al. 2005 

) and 

possi- 

ble planetary companionship (IVlullally et a 

T Tt 11/^ /~\T Tr\lllV^/Y TTTl T" Ll 1 T" Lc Ia If’ 

. 2008 
■jkrTjni 

)• By 

O 


become the most extensively observed cool-edge ZZ Ceti. 
That this outburst behavior has not been previously de¬ 
tected in other cool-edge pulsators may be due partially 
to this selection bias. The routine practice of dividing 
low-order polynomials from ground-based light curves to 
correct for changes in atmospheric extinction may also 
have obscured this behavior in previous observations. 


4. WHITE DWARF PULSATION ANALYSIS 

The 20-month light curve yields one of the richest pul¬ 
sation spectra ever resolved for a ZZ Ceti variable. If 
the energetic outbursts do, in fact, originate from the 
white dwarf, they likely affect the stellar pulsations in a 
measurable way. We explore the observed variations in 
pulsation mode properties that may be related to these 
outbursts, but unfortunately the signal-to-noise of the 
Kepler light curve is insufficient to study these variations 
on timescales less than the mean outburst timescale in 
much detail. We proceed with an asteroseismic analysis 
that assumes that the outbursts do not significantly alter 
the locations of the peaks in the power spectrum of KIC 
4552982 and reserve comments on the possible interplay 
between these dynamic processes for our concluding re¬ 
marks. 


4.1. Mode Stability and Frequeney Determination 

Figure!^ shows the Fourier transform (FT) of the entire 
Qll-QlTTight curve through the full region of significant 
pulsational variabilit y. This FT was comp uted with the 
Period04 software (Lenz & Breger 2004). We exclude 
all known instrumental artifacts that are harmonics of 


the long-cad ence sampling rate of 566.41 /iHz (Gilliland 
et al.||2010b|) from this analysis. The power at the low- 
frequency limit of the FT is introduced primarily by the 
aperiodic outbursts and mostly goes away if we remove 
the outburst events from the light curve before comput¬ 
ing the FT. Since this 1// noise decays sufficiently before 
the low end of the frequency range of pulsations and the 
FT of the full light curve (including outbursts) yields an 
overall lower noise level in the spectrum, we do not ex¬ 
clude the outbursts from the light curve for our Fourier 
analysis. 

The features of the FT of particular asteroseismic sig¬ 
nificance are the 17 bands of power in the 600-1450 
yuHz frequency range and the sharp rotational triplet fea¬ 
ture surrounding 2765.66 /iHz. The broad power bands 
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Frequency (//Hz) 

Fig. 8.— Fourier transform of the entire Q11-Q17 Kepler light curve of KIC 4552982 in the region of significant pulsational variability 
(amplitude scale in mma = 0.1 %). 


0 


1 


2 


3 


4 


5 


6 


7 


8 Amplitude [mma] 


56400 

56350 

56300 

56250 

_ 56200 

o 

o 

o 

i 56150 

csj 

Q 

E 56100 

O) 

E 

i- 

56050 

56000 

55950 

55900 

55850 



700 800 900 1000 1100 1200 1300 2740 2760 2780 

Frequency [jiHz] 


Fig. 9. — Sliding FT of a 5-day window over the entire Kepler light curve. The color bar at the top of the figure gives the amplitude 
scale in mma (= 0.1 %). The modes in the 700-1300 yuHz frequency range (left panel) are observed to wander in phase (frequency) and 
amplitude while the mode at 2757.54 iiHz is relatively stable. 


clearly do not reflect the sharp, 0.020-//Hz-wide spectral 
window (the signature of a perfect sine wave sampled 
identically to the data). This is the result of modulation 
of the pulsation properties of the observed modes. This 
modulation in amplitude and frequency is clearly seen in 
the sliding FT of Figure The sliding FT is computed 
by sliding a 5-day window over the entire light curve and 
computing the FT at 1-day steps. This time-resolved 2- 
dimensional FT is shown in color contour, where darker 
areas indicate greater power in the FT. We observe that 
all the pulsations in the 700-1300 //Hz range (left panel) 
show rapid changes in frequency and amplitude on a 
timescale of a few days. The frequency shifts of the main 


peaks appear to be correlated, demonstrating that these 
drifts cannot be ascribed purely to stochastic behavior of 
the individual modes. The higher amplitudes in the slid¬ 
ing FT compared to the full FT in Figure [^demonstrates 
that the instantaneous mode amplitudes are larger than 
the peaks of the widened power bands. These data of¬ 
fer the most extensive coverage of m ode variations that 
are typical of cool ZZ Cetis (see, e.g., Pfeiffer et al.|1996 


Kleinman et al. 1998 Hermes et al.|2014]). 'The variations 

lie o 


in mode properties may be related to the outbursts ob¬ 
served in the light curve (Section]^, but we were unable 
to demonstrate this connection from the data. Mean¬ 
while, the mode at 2757.54 //Hz (right panel) is relatively 
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Fig. 10. — Top: Lorentzian fits to 17 bands of pulsational power in the 4-/i,Hz boxcar smoothed FT. The horizontal line at 0.0468 mma 
is our significance criterion. The Lorentzians are best fits to the power (amplitude squared) above the median noise level in hand-selected 
(black) regions of the light curve. Bottom: We refine our fit parameters by simultaneously fitting these 17 bands to the original, unsmoothed 
FT. The final parameters of these Lorentzians are listed in Table 


stable. 

We interpret each band of power to be linked to a sin¬ 
gle significant pulsation mode in the star. Since the wide 
nature of these bands prevents us from simply selecting 
and prewhitening the highest peaks in the FT, we in¬ 
stead adopt a method of fitting each band in the power 
spectrum (the squared FT) with a Lorentzian function: 


/(0 = 


p^2 


{v - VoY + 72 


+ {FT{vf). 


( 1 ) 


While Lorentzians are used to fit signatures of the 
stochastically driven pulsations in Sun-like stars and red 
giants, we emphasize that stochastic driving is unlikely to 
be efficient in pulsating white dwarfs since the stochastic 
driving timescale 1 second) is much shor ter than th e 
observed pulsation periods 10 minutes; Saio 2013). 
We choose to fit Lorentzians as a convenience! We 
use slightly different methods for identifying and fitting 
Lorentzians to the significant frequencies in each region 
of the FT that are adapted to suit the different feature 
widths. 

For the dense bands of pulsational power in the 600- 
1450 /iHz range, we first identify which bands are signif¬ 
icant from the 4-//Hz boxcar-smoothed FT displayed in 
the top panel of Figure pT| We use a bootstrap method to 
calculate a very conservative significance threshold. We 
randomly rearrange the light curve points 10,000 times. 


keeping the same time sampling. We then calculate the 
FT in the 600-1450 //Hz range, apply the 4 //Hz smooth¬ 
ing, and record the highest value in the smoothed FT 
for each shuffled light curve. We set our significance 
threshold equal to the 99.7 percentile in the distribu¬ 
tion of maximum values in each shuffled, smoothed FT. 
We find this value at 0.0468 mma (1 mma = 0.1%) in 
the smoothed FTs. This means there is only a ^ 0.3% 
chance that we would find any peak in the smoothed 
FT above 0.0468 mma due to noise alone. However, in 
regions near a pulsational power band, the wide feature 
can raise actual noise peaks above this threshold, so we 
conservatively fit only unambiguously significant features 
near larger power bands. The peaks of the features we 
select are each marked with a x in the top panel of Fig¬ 
ure!^ Our least squares Lorentzian fits use 5 //Hz ini¬ 
tial guesses for the HWHM at these locations (selected 
so that the resulting fits match the intended features by 
eye). We fit Lorentzians only to the regions immediately 
surrounding the selected features (the black portions of 
the smoothed FT in the top panel of Figure 10). The 
resulting Lorentzian fits to these features are olsplayed 
in the figure. 

While the smoothed FT is useful for selecting signifi¬ 
cant peaks in frequency space, the smoothing can affect 
the HWHMs of the Lorentzian fits. Since these HWHMs 
may have an astrophysical significance, we refine our fit 
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TABLE 2 

Properties of Lorentzian fits to significant 

FREQUENCIES OF PULSATIONAL VARIABILITY IN THE FT. 


Mode 

Period 

(s) 

Frequency 
{uq; /iHz) 

HWHM 
( 7 ; mHz) 

Lorentzian height 
(P; mma^) 


1498.32 

667.42 

0.34 

0.0062 

h 

1362.95 

733.70 

0.26 

0.0056 

/3" 

1333.18 

750.09 

0.46 

0.0050 


1301.73 

768.21 

0.97 

0.0070 


1289.21 

775.67 

1.38 

0.0132 


1244.73 

803.38 

4.50 

0.0023 

fr^ 

1200.18 

833.21 

2.15 

0.0018 


1158.20 

863.41 

6.59 

0.0055 

h 

1100.87 

908.38 

4.51 

0.0023 

/lO 

1053.68 

949.06 

8.61 

0.0032 

/ii 

1014.24 

985.96 

7.13 

0.0066 

/l2 

982.23 

1018.09 

2.21 

0.0081 

/l3" 

950.45 

1052.13 

4.58 

0.0246 

/l4" 

907.59 

1101.82 

4.75 

0.0189 


866.11 

1154.59 

3.68 

0.0266 

fl6‘^ 

828.29 

1207.31 

4.67 

0.0202 

/17" 

788.24 

1268.65 

3.93 

0.0029 


362.64 

2757.54 

0.013 

0.8743 


361.58 

2765.66 

0.007 

0.0260 


360.53 

2773.71 

0.009 

0.0260^ 


^ Likely part of the i = \ sequence based on alignment with the 
mean 1=1 period spacing or observed rotational splitting. 

^ Likely rotationally split components of an ^ = 1 mode. 

This fit matches the single peak for this component in the 
power spectrum with a height corresponding to peak power. 


parameters by fitting to the unsmoothed FT. We refit all 
of the Lorentzians concurrently with a common bias-level 
parameter to fit the average noise level. The result of 
this fit is plotted over the unsmoothed FT in the bottom 
panel of Figure Because these power bands encom¬ 
pass a forest of high- and low- amplitude peaks, the final 
Lorentzians fit through the average amplitudes and do 
not reach the maximum peaks of the features. The cen¬ 
tral frequencies, HWHMs, and heights of the Lorentzian 
fits are listed with the central mode periods in Table 
Their relation to Equation His also indicated. 

Figure shows a higher-mequency region of the origi¬ 
nal, unsmoothed FT. Here we observe a triplet of signif- 



FlC. 11.— The rotational triplet surrounding the 2765.66 p,Hz 
mode. The power was fit with the displayed set of Lorentzians. 
The highest-frequency component consists of a single point above 
the noise that is fit to its peak power by the Lorentzian. 


leant variability with near-even frequency spacing. We 
calculate a significance threshold of 99.7% confidence at 
0.167 mma by a bootstrap approach similar to our calcu¬ 
lation for the lower-frequency power bands. The differ¬ 
ence between the 0.167 mma significance threshold here 
and the 0.0468 mma threshold for the smoothed FT at 
lower frequency results purely from the boxcar smooth¬ 
ing of the latter. While only the two lowest-frequency 
components rise above this significance threshold, we re¬ 
lax our criterion for the third peak because it falls in step 
with th e oth ers as is expected for rotational splitting (see 
Section 4.4) and is only barely below our very conser¬ 


vative sigmfiance threshold. We simultaneously fit the 
three features with Lorentzians as we did for the wider 
power bands using 1 /iHz for the initial guesses for the 
HWHM and adopting the bias level fit in the previous 
region. The third component consists of a single peak 
above the noise, and the height of the Lorentzian fits the 
peak power of that mode. These best-fit parameters are 
included in Tabled 

We adopt these 20 frequencies as the mean pulsation 
mode frequencies of KIC 4552982. Since we used an un¬ 
conventional method for determining these mode param¬ 
eters for pulsating white dwarfs, it is difficult to assign 
absolute uncertainties. We advise that researchers wish¬ 
ing to fit their own asteroseismic models to these fre¬ 
quencies weight their fits by 1/HWHM^ for each mode. 
Since we resolve the frequency variations over our ex¬ 
tended light curve, our frequency determinations should 
be far more accurate than the HWHM. 

4.2. Period Spacing 

The periods of the nonradial ^f-mode pulsations char¬ 
acteristic of white dwarfs of a given degree {£) and se¬ 
quential radial order {k) reach an even spacing at the 
asymptotic limit of high k. Because of geometric cancel¬ 
lation, we exp ect the pulsations w e detect to be of degree 
^ = 1 or ^ = 2 (Dziembowski 1977). With such a rich pul¬ 
sation spectrum, we can hope to sample these sequences 
sufficiently well to determine the mean period spacing 
with radial order of one or both of these sequences. 

We approach this measurement by taking the period 
transform of the period transform of the entire light curve 
as is described in its application to the D OV (pulsating 
hot pr e-white dwarf) star PG 1159-035 in Winget et al. 
(1991, Section 5.2.1). We emphasize that this approa ch is 
independent of the period determination in Section |4.1[ 
We arrive at the period transform by simply inverting 
the x-axis of the FT in the region of dense pulsational 
power bands (between 500 and 1500 /iHz, corresponding 
to a period range 667-2000 s). This intermediate result 
is shown in Figure If the signals from pulsations 
correspond to one or more sequences with evenly spaced 
periods, the period transform of this period transform 
should reveal the spacings. We show the resulting power 
spectrum in Figure and mark two peaks of interest: 
a significant spacing at 41.9T0.2 s and its overtone at 
20.97 ± 0.02 s (uncertainties determined from 100 Monte 
Carlo simulations with Period04). This overtone arises 
from the decidedly non-sinusoidal nature of the periodic 
signal in the period transform. Using the phase informa¬ 
tion we get from least-squares sinusoidal fits, we mark the 
locations of expected modes in the even period-spacing 
sequence as vertical dotted lines in Figure [l^ and indicate 
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Period (s) 



Mean Period Spaeing (s) 


Fig. 12.— The period transform of the full Kepler light curve 
in the region of observed pulsational power bands. The locations 
of the expected i = 1 modes from the mean period spacing deter¬ 
mined in Figureare marked as dotted vertical lines. These lines 
are drawn darker where they fall within 7 seconds of one of our 
measured mode periods listed in Table 



Fig. 14.— Histograms of mean period spacings calculated for 
i = 1 and i = 2 sequences throughout our model grid. The 
measured spacing for KIC 4552982 of 41.9 seconds fits exclusively 
within the i = 1 distribution. 



Teff (K) 


Fig. 13.— The period transform of the period transform reveals 
a significant mean period spacing at 41.9 seconds and its harmonic. 
This likely corresponds to a mean period spacing oi I = 1 modes. 

those modes that fall within 7 seconds of the measured 
spacing with darker dotted lines and with footnotes in 
Table |2l 


Fig. 15.— The locations of models through M* - Tgff space with 
mean 1=1 period spacings equal to the measured spacing of 41.9 
seconds. The individual runs of solutions correspond to log Mh/M* 
values from -4.00 (bottom track) to -6.00 (highest track) in 0.50 
resolution. All models that fit the observations have M* > 0.60 
Mq . The spectroscopically determined parameters are overplotted 
in white, and are in slightly better agreement with models nearer 
the high-Mn end. 


4.3. Comparison with Asteroseismic Models 

We consider a grid of more than 14,000 DA white dwarf 
cooling models that we c alculated from the White Dwarf 


Evolu tion Code (WDEC;|Lamb fc van Hor n|1975||Wood 
1990|) followi ng the treatment described in Bischoit-Kim 


et al. (2008). We vary the following three parameters 


that most influence the mean period spacing, with noted 
resolution: 10,000 K < Teff < 12,000 K (200 K reso¬ 
lution), 0.500 Mq < < 0.800 Mq (0.005 Mq reso¬ 

lution), and —6.00 < logMn/M* < —400 (in steps of 
0.10). This set of models widely encompasses the spec¬ 
troscopic M^, Teff values of KIC 4552982. Each of these 
models has a helium layer mass of logMne/M^ = —2.00 
(cl ose to the value calculated f or 0.7 Mq white dwarfs 
by |Lawlor & MacDonald 2006). The core profiles have 
central abundances of 307o carbon and 70% oxygen and 


are homogenous out to 0.5 M^/M^ (chosen to be in good 
agreement with S alar is et al. 1997|). 

The measured mean period spacing of 41.9 ± 0.2 s fits 
well within the distribution of calculated mean spacings 
in our models for ^ = 1 modes. Eor i = 2 modes, the 
models predict mean spacings nearer to 26 s. These dis¬ 
tributions are compared in Eignre[T^ . 

We suggest that the modes thatTall in step with the 
measured mean period spacing are likely part of the i = 1 
sequence with footnotes in Table Undoubtedly, some 
£ = 2 modes may by chance fall in line with this se¬ 
quence and be misidentified. Mode trapping and the 
dominance of different rot at ionally split components of 
the modes (see Section 4.4) can shift individual frequen¬ 
cies significantly from an equal period spacing, so some 
of the modes that we do not observe to closely match the 
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mean spacing may also have ^ = 1. The dominance of 
different rotationally split compone nts should not p eatly 
affect our asteroseismic inferences ( Metcalfe]|2QQ3| . 

We display the run of models with mean = 1 spac- 
ings of 41.9 s through M^, Teff space for selected hydro¬ 
gen layer masses in our model grid in Figure The 
mean period spacing of ^-modes is most sensitive to the 
overall mass, and our model comparison demonstrates 
asteroseimically that KIC 4552982 is likely more massive 
than 0.6 Mq — in good agreement with the spectroscopic 
mass determination. 

The spectroscopic point is plotted with error bars 
over the logMn/M^ tracks for reference. We observe 
that while the uncertainty in the spectroscopically de¬ 
termined spans multiple tracks at the resolution of 
our model grid, the spectroscopic point agrees best with 
a logMH/TA value of -4.70 and generally falls nearer 
to the maximal logMn/M^ extent of our model grid at 
-4.00. 

4.4. Rotational Splitting 

At higher frequency than the region of dense pulsa- 
tional power bands, we noted three modes split evenly in 
frequency (FigurepTj). Even frequency spacing can result 
from r otational splitting of a mode into 2^+1 compo¬ 
nents (Unno et al.||1989 ). While we emphasize that our 
absolute uncertai ntie s m the individual frequencies deter¬ 
mined in Section [4T] are smaller than the HWHM of our 
Lorentzian fits, we will treat the HWHM here as proxies 
for uncertainty. From these values, we calculate a mean 
frequency splitting of 8.09 ± 0.02 //Hz. 

With only three components in the rotational multi- 
plet, we are unable to definitively assign a spherical de¬ 
gree, to this mode, so we consider the three most likely 
possibilities: these are either all m = — 1 , 0,1 compo¬ 
nents of a ^ = 1 triplet, three consecutive components of 
a ^ = 2 quintuplet, or the m = — 2 , 0 , 2 components of a 
£ = 2 quintuplet. 

We derive first-order analytical rotation rates under 
the assumption that the rotation period is much longer 
than the period of the split mode and that the white 
dwarf rotates as a solid body. We use the rela tion that 
^^k£m = i'^/Prot) X (1 “ Cm) (|Brickhill||1975|) for solid 
body rotation. At the limit of high k [k ^ lU), Cm ap¬ 
proaches l/^(-^ + 1). If we assume that this is an = 2 
mode, we are safely in the asymptotic limit {k ^ 11 or 
12 according to our models). If these are adjacent split¬ 
tings of a ^ = 2 mode (i.e., m = —2, —1, 0, m = —1, 0,1 
or m = 0,1, 2), we derive a rotation period of 28.63±0.07 
hours. If these are £ = 2, m = —2,0,2 modes, we get 
twice that rotation period at Prot = 57.26 ± 0.13 hours. 
If this is the rotationally split triplet of a ^ = 1 mode, 
we cannot safely assume that these modes are in the 
asymptotic limit. We guide our interpretation for this 
case by comparing the mode periods to our asteroseismic 
model with parameters Teff = 10,800 K and = 0.700 
Mq and log Mu /= —4.70. This model most closely 
matches our measured mean £ = 1 period spacing and 
our spectroscopically determined atmospheric parame¬ 
ters. The measured 361.58-s period measured for the 
central component of the triplet falls between the follow¬ 
ing two £ = 1 modes in the model: /c = 5 at 311.6 s with 
Cm = 0.491, and /c = 6 at 379.0 s with Cm = 0.492. 
These Cm values essentially match the asymptotic re¬ 


sult, and we adopt the mean of these values to calculate 
a rotational period of 17.47 ± 0.04 hours. These periods 
are all plausible considering the range of asteroseismi- 
cally determine d rotation periods of other white dwarfs 
(Kawaler 2004), and they support our assumption that 
the rotation period is much longer than the pulsation 
periods. 

Revisiting the region of wide puls ational power bands 
between 600-1450 /rHz (Figure ]1Q|), we see suggestions 
of structure at a similar ^ 8 /iHz spacing. This includes 
a pair of frequencies that passed our significance crite¬ 
ria and are indicated with a footnote in Table H /4,/5. 
Since these bands fall in line with the £ = 1 m ean pe- 
riod spacing pattern from our analysis in Section [T^ we 
prefer the £ = 1 rotation period of 17.47 ± 0.04 hours for 
this star. 

5. CONCLUSIONS AND FUTURE WORK 

The Kepler light curve for KIC 4552982 exhibits two 
features of great interest: clear outburst phenomena and 
a rich spectrum of pulsations with frequency and ampli¬ 
tude modulations. We argue that these outbursts likely 
originate from the pulsating white dwarf and are the first 
observations of a new astrophysical phenomenon. Such 
energetic events likely affect the pulsations in a measur¬ 
able way. Owing to the low signal-to-noise of the Kepler 
photometry of this faint target, we are unable to study 
changes in the pulsation spectrum strictly before and af¬ 
ter the observed outburst events. 

The most compelling evidence for the outbursts affect¬ 
ing the pulsations is the relative sharpness of the triplet 
of high-frequency modes surrounding 2765.66 /iHz com¬ 
pared to the wide bands of power at lower frequency. 
If these modes are all of the same degree, £, the lower- 
frequency, higher-/c modes are more sensitive to regions 
nearer the surface of the star. We demonstrate that most 
of the modes we observe are likely part of a ^ = 1 se¬ 
quence. The detected relative unsteadiness of the lower- 
frequency modes may then be an indication that the out¬ 
bursts are a surface phenomeno n (a similar a rgument is 
made for solar-type pulsators by Karoff||2Q14 ). However, 
variations in the pulsation frequencies and amplitudes 
are commonly observed in cool ZZ Cetis, so this is not 
the only viable interpretation. If outbursts of the type 
we see in KIC 4552982 turn out to be common of cool ZZ 
Cetis, they may be related to the frequency modulations 
observed in these other stars. 

We updat ed the spectroscopic parameters of |Hermes| 
et al. (2011) with corrections from 3D convective simu- 
lations to get Teff = 10,860 ± 120 K, log^ = 8.16 ± 0.06, 
and = 0.69 ± 0.04 Mq. These parameters place KIC 
4552982 at the extreme cool edge of the ZZ Ceti insta¬ 
bility strip where pulsations are just shutting down. We 
demonstrate that the average duration (« 10 hours) and 
recurrence (~ 2.7 days) timescales observed for the out¬ 
bursts are of order the thermal timescale at the base 
of the convection zone and suggest that these outbursts 
could be a feature common to all white dwarfs with this 
convection zone depth, and therefore temperature. 

We identify 20 independent pulsation frequencies in¬ 
cluding rotationally split components as characterized in 
Table Comparing the measured mean period spacing 
of the nearly sequential /c-overtones of likely I = 1 modes 
to our models, the asteroseismology supports the spec- 
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troscopic finding for the white dwarf mass: M* > 0.6 
Mq. These spectroscopic parameters match our astero- 
seismic models with 41.9 i 0.2 s mean period spacings 
slightly better closer to the thick outer hydrogen layer 
end of our grid (logMn/M^ = —4.00) rather than the 
thin limit of our models (logMn/M^ = —6.00). We also 
derive a likely rotation period of 17.47 ±0.04 hours if the 
observed modes surrounding 2765.66 /iHz are a rotation- 
ally split £ = 1 triplet. 

The sheer extent of this nearly continuous 20-month 
light curve has provided an unparalleled look at dy¬ 
namic white dwarf processes. Continued asteroseismic 
interpretation of the identified pulsation frequencies and 
outburst characteristics will hopefully lead to additional 
constraints on the internal structure of KIC 4552982. 
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